NUWC  NPT  Technical  Document  10,006 
3  March  1992 


Ternary  Weak-Signal  Detection 
in  Non-Gaussian  Noise: 

A  Preliminary  Analysis  for 
"Ho:  N  vs  Hi:  N  +  Si  vs  Hr.  N  +  S2" 
With  Independent  Sampling 

D.  Middleton 
Sonalysts  Inc. 


J.  G.  Kelly 

Weapon  Systems  Technology  and  Assessment  Department 


DTIC 


ELECTE 
MAY  2  4 1993 


E 


Naval  Undersea  Warfare  Center  Division 

Newport,  Rhode  Island 


£3  5  20  06 1 


93-11510 

_ . 


Approved  for  public  release;  distribution  is  uninled. 


-  A 


PREFACE 


This  document  was  prepared  for  the  Weapon  Systems  Technology  and 
Assessment  Department  under  Contract  N66604-88-D-0090-0021. 


Reviewed  and  Approved:  3  March  1992 


I^LUaJaaM 

F.  L.  White 

Head,  Weapon  Systems  Technology 
and  Assessment  Department 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-0188 


tuHtoi  fur  nArihii  M  iai  fa  hi 


.  todudng  np—ww  PanUudnp  Ma 
HMiMpMVV*  BMMNl,aPI«*iOtaP 


i  m  any  i 

,  mraaana*  far  MonaaSan  < 

jcPan  Puna  yw^PWfc  WaaN wptaw.  PC  : 


1.  AGENCY  USE  ONLY  (Isaim  Blank) 

2.  REPORT  DATE 

3  March  1992 

J.  REPORT  TYPE  ANO  OATES  COVERED 

Final 

4.  TITLE  AND  SUBTITLE 

Ternary  Weak-Signal  Detection  in  Non-Gaussian  Noise:  A  Preliminary  Analysis  for 
"Hq:  NvsHj:N  +  Sjvs  H2:  N  +  S2"  With  Independent  Sampling 

S.  FUNDING  NUMBERS 

6.  AUTHORS) 

D.  Middleton 

J.  G.  Kelly 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADORESS(ES) 

Naval  Undersea  Warfare  Center  Division 

Newport,  Rhode  Island  02841-5047 

A  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

TD  10,006 

9.  SPONSORINGMONfTORING  AGENCY  NAME(S)  AND  ADORESS(ES) 

10.  SPONSOR!  NQ/WONnrOFUNa 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 


12a.  OISTRIBUTIOWAVAILABIUTY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (tUximum  200  words) 

A  general  analysis  of  the  Ternary  Gass  (M  =  2):  H  0:  N  vs  H  j:  S  j  +  N  vs  H2:  S2  +  N  of  signal  detection  problems  is 

presented,  for  completely  general  signals,  i.e.,  both  broad-  and  narrow-band,  deterministic  or  random,  in  generalized  (i.e., 
non-Gaussian)  noise,  in  the  limiting  threshold  regime.  This  includes  optimum  threshold  algorithms  and  system  performance, 
as  measured  by  the  appropriate  error  and  detection  probabilities.  The  present  treatment,  however,  is  subject  to  the  following 
constraints:  (1)  independent  noise  sampling;  (2)  ambient  noise  models,  i.e.,  noise  independent  of  the  signals;  (3)  uniform  cost 
functions,  e.g.,  C0  (>  0)  for  errors,  and  C  j  =  0  for  correct  decisions.  Under  these  conditions,  only  three  principal  parameters 

are  needed:  ol2,  o22  =  signal  detection  parameters  (=  "output  (S/N)  2")  and  the  correlation  coefficient  p12  (=  p)  between  the 
two  (threshold)  test  statistics  (or  detection  "algorithms")  Z  j,  Z2,  apart  from  the  a  priori  probabilities  (q,  p  j,  p2)  of  the 
presence  of  noise  alone,  S  j,  and  S2.  Next  steps,  toextend  the  treatment  to  the  general  case  (M  =  3):  H  j:  N  +  Sj,  vs 


H2:  S2  +  N  vs  H3:  S3  +  N,  and  to  include  correlated  noise  samples,  are  noted. 


14.  SUBJECT  TERMS 

Ternary  Detection  Coherent  and  Incoherent  Reception  General  Signals 

Threshold  Signal  Detection  Generalized  Noise 

IS.  NUM8ER  OF  PAGES  5Q 

IS.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 

1«.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

SAR 

NSN  7540-0 1-280-S600 


Standard  Form  2M  (Rav.  2-S8) 
fiaacHbad  by  ANM  W  MMI 
in-tn 


TABLE  OF  CONTENTS 


Section  Page 

LIST  OF  ILLUSTRATIONS .  iii 

LIST  OF  PRINCIPAL  SYMBOLS .  iii 

1  INTRODUCTION .  1-1 

2  TERNARY  DETECTION:  Ho:  NvsHi:  N+S}  vsH2:  N+S2:  GENERAL 

THEORY,  WITH  EQUAL  COST  ASSIGNMENTS .  2-1 

2.1  The  Decision  Process,  (M=2) .  2-1 

2.2  BayesRisk .  2-3 

2.3  Canonical  Threshold  Detection  Algorithms,  Si  and  S2 .  2-5 

2.3.1  Coherent  Detection  <s>  *  0:  (k  -  1,2) .  2-5 

2.3.2  Incoherent  Detection  <s>  *  0:  (k  =  1,2) .  2-6 

2.4  Performance  Probabilities:  Ternary  (M  =  2) .  2-8 

3  EVALUATION  OF  THE  PERFORMANCE  PROBABILITIES: 

TERNARY  (M  =  2) .  3-1 

3.1  The  Normalized  PDFs  of  Z .  3-1 

3.1.1  Coherent  Detection . 3-2 

3.1.2  Incoherent  Detection .  3-2 

3.2  The  Performance  Probabilities:  (M  =  2) .  3-4 

3.2.1  Probabilities  of  Correct  Detection .  3-4 

3.2.2  Probabilities  of  Error .  3-9 

4  ADDITIONAL  SPECIAL  CASES .  4-1 

4. 1  (<H,2  — >  <r2fi  <«») .  4-1 

4.2  The  Cases  p  =  ±  1 .  4-2 

4.3  TheCasep  =  0 .  4-4 

4.4  Reduction  to  the  Binary  Cases:  Hi  vs  H2  and  Hi  vs  Hq .  4-5 


5  CONCLUSIONS .  5-1 _ 

APPENDIX  A  -  EVALUATION  OF  TERNARY  PROBABILITY  rf 

INTEGRALS .  A-l  □ 


APPENDIX  B  -  CORRECTION  FOR  SOME  RESULTS  IN 
REFERENCE  3 .  B-l  — 


REFERENCES .  R-l  odes 


I  I  Avn:! 

DTIC  QiiaLrfV  5 


i/ii 

Reverse  Blank 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

2-1  Decision  Regions:  Yi,  Y2  and  Boundaries .  2-2 


LIST  OF  PRINCIPAL  SYMBOLS 


A)c  =  derector  bias 

a(^  =  (conditional)  error  probs,  k  =  1,2: 
"noise  alone"  =  true  state 
=  (normalized)  signal  amplitudes 

PT* etc- =  (conditional)  error 
probs. 

C0  =  cost  of  errors 


/j  -  equation  (2- 1 2c)  nonlinear 

functions 

A  1,2  =  likelihood  ratios 
)  =  costs 

M  =  no.  of  signal  classes;  total  no.  of 
spatial  samples 

Pk  =  ratio  of  a  priori  probs.,  k  =  1,2 
m  =  111th  spatial  sample  (sensor) 


e,  £q  =  signal  epochs 

gOO*  =  opt  threshold  algorithms  (k  = 

1,2) 

YO.1,2  -  decisions 

Ho,  Hi,  H2  =  Hypothesis  states 

Ic  =  incomplete  T-function 
1(2)  =  Basic  2nd-order  probability 
integral,  equation  (A.l-1). 

J  =  MN  =  total  no.  space  x  time  samples 
J( )  =  equation  (A.  l-5a)  et  seq.,  a 
definite  integral 

K.(i),(2)  =  thresholds 


N  =  noise;  total  no.  of  time  samples 
n  =  n*h  time  sample 

Ci\j  =  Doppler  (angular)  frequency 

P(q)  =  probs.  of  correct  detection 
p^)  =  (conditional)  prob.  of  correct 
detection 

V  =  total  background  noise  intensity, 
equation  (2- 12c) 

(j><m)  =  order  derivative  of  error 
function 

Q2,  q2  =  2nd-order  PDFs,  noise  alone 

R^,  =  Bayes  risk 

p,  pl2  =  correlation  coeff.  for  Z\Z  1 


L(2),L<4>  =  noise  statistics,  equations 
(2- 12c),  (2- 13b) 


LIST  OF  PRINCIPAL  SYMBOLS  (Cont  d) 


s®=  (normalized)  signal  sample  W2  =  PDF  of  random  variables  (xi,X2) 

o£  -s  detection  parameter,  k  =  1,2 

Zj,  Z2  =  data-dependent  part  of  detector, 
ATm=  array  time  delay,  equation  (2-  14b)  see  equation  (2-12) 

0(j  \  6,^  =  (normalized)  signals,  z,  21,2  =  normalized  random  variables, 

k  =  1,2  equation  (3-1) 

0  =  error  function 


tv 


TERNARY  WEAK-SIGNAL  DETECTION  IN  NON-GAUSSIAN  NOISE: 

A  PRELIMINARY  ANALYSIS  FOR 

"Hq:  N  vs  Hi:  N  +  Si  vs  H2:  N  +  S2"  WITH  INDEPENDENT  SAMPLING 


1.  INTRODUCTION 

An  important  extension  of  the  usual  binary  detection  theory  involving  generalized, 
or  non-Gaussian  noise  (see  reference  1,  for  example)  is  the  ternary  case  involving  two 
signals  (M  =  2)  in  noise,  as  well  as  noise  alone.  The  basic  hypothesis  states  are  Ho:  N, 
"noise  alone,”  vs  Hj:  N  +  Si,  "signal  1  in  noise"  vs  H2:  N  +  S2,  "signal  2  in  noise." 

This  is  a  special  case  of  the  general  M-ary  detection  model,  described  in  section  23.1  and 
outlined  in  (3),  (b)  of  section  23.1  of  reference  2  for  the  cases  M  =  2,3.  (Also  see 
reference  3  for  some  recent  ternary  detection  analysis  involving  purely  Gaussian  noise, 
where  effective  erasure  of  magnetic  tapes  is  considered.) 

The  present  study  provides  new  results:  (1)  a  canonical  formulation  for  these  Hq  vs 
Hi  vs  H2  cases,  for  arbitrary  signals,  e.g.,  broadband  as  well  as  narrowband,  deterministic 
or  purely  random,  with  (2)  spatial  as  well  as  temporal  sampling — thus  including  the 
potential  for  adaptive  beam  forming  or  generalized  "matched  field  processing"1  here,  as 
well  as  the  usual  purely  temporal  cases,  involving  preformed  beams;  and  (3)  threshold 
reception  in  generalized  noise,  which  in  turn  may  be  fully  described  by  Class  A  and  B 
models  in  most  instances.  (See  references  4  through  6.  For  a  review  of  these  latter  in 
related  telecommunications  applications,  see  references  7  and  8.) 

The  limitations  of  the  present  analysis  are  twofold:  (1)  restrictions  to  independent 
noise  sampling,  and  (2)  the  simple,  but  reasonable  cost  assignments  that  provide  the  same 
costs  to  "failure,"  i.e.,  wrong  decision,  and  identical  (but  lesser  costs,  necessarily)  to 
"success,"  or  correct  decisions.  The  latter  is  usually  not  a  serious  restriction  in  most  cases, 
and  it  greatly  simplifies  the  analytic  specification  of  the  decision  regions  (see  equation 
(23.20)  and  figure  23.3  of  reference  2  vis-i-vis  the  case  (b)  in  section  23.1  cited  above). 
The  former  simplification,  however,  can  lead  to  noticeably  suboptimum  performance  vis-4- 
vis  correlated  noise  sampling,  when  the  latter  can  or  should  be  applied.  (This  limitation 
will  be  removed  in  a  subsequent  study,  along  the  lines  of  references  9  and  10.) 


1  This  term  is  used  for  the  spatial  analogue  of  the  earlier  concept  of  "matched  filter  processing"  for  the 
temporal  representations  of  the  signal  (output  of  receiver  array  or  preformed  beam).  This  is  a  local  concept, 
and  not  a  global  one,  in  the  sense  of  the  more  recent  "Matched  Reid  Processing"  (MFP)  or  passive  signal 
locations  (PSL)  developed  elsewhere  (see  reference  17,  for  example). 
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The  generic  application  of  the  present  analysis  is  to  those  situations  in  signal 
detection  where  one  of  the  signals,  i.e..  Si,  is  the  "desired"  signal,  and  S2  embodies  an 
unwanted  signal.  For  the  current  application  to  (active  or  passive)  underwater  acoustic 
systems,  Sj  can  represent  a  target  return  for  an  emitting  source  whose  presence  must  be 
determined,  while  S2  is  a  possible,  confusing  or  un desired  ambient  emission,  a  false  target 
or  decoy.  At  the  same  time,  the  non-Gaussian  character  of  the  background  noise  arises 
from  a  variety  of  physical  mechanisms:  (1)  isolated,  undesired  signal  sources,  e.g., 
shipping,  biological  emitters,  other  man-made  and  natural  phenomena,  etc.;  and  (2) 
distributed  sources,  such  as  reverberation,  ice  noise,  surface  effects,  e.g.,  waves,  rain,  etc. 
Reverberation  noise,  or  scattering,  is  by  definition  signal  dependent,  as  opposed  to 
"ambient  noise,"  whose  sources  are  outside  the  control  of  the  transmitter  in  question  (see 
reference  1 1).  Ignore  for  the  present  the  signal-dependent  nature  of  reverberation, 
remembering,  nevertheless,  that  its  detailed  statistical  structure  depends  on  the  original, 
injected  signal,  as  well  as  on  the  particular  scattering  mechanism,  e.g.,  ocean-air  interfaces, 
ocean  bottoms,  etc.  In  any  case,  regardless  of  the  mechanism  involved,  the  canonical 
models  here  are  applicable  not  only  in  the  present  underwater  acoustical  context,  but  in  the 
more  usual,  analogous  telecommunications  case,  as  well  (see  reference  12). 

This  paper  is  organized  as  follows:  section  2  provides  (1)  a  general  formulation  of 
the  ternary  "two-signal"  (M=2)  detection  cases,  for  the  simple  but  common  cost 
assignments  Cq  (>  0):  "incorrect",  and  C  =  0:  "correct;"  and  (2)  specialization  to 
threshold  signals  and  independently  sampled  noise  fields.  Section  3  gives  the  evaluations 
of  the  associated  performance  probabilities  in  the  threshold  regime.  Various  limiting  and 
special  cases  are  discussed  in  section  4,  which  also  serve  to  verify  the  general  results  of 
section  3.  The  paper  concludes  in  section  5  with  a  short  summary  of  the  scope  and 
limitations  of  the  present  analytical  model,  and  includes  a  series  of  next  steps  to  be  pursued 
in  the  mathematical  and  numerical  analyses  of  the  general  ternary  case.  Numerical  results, 
requiring  alternative  series  development,  are  reserved  for  a  subsequent  report  (see  reference 
13). 
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2.  TERNARY  DETECTION:  HO:  N  vs  Ht:  N  +  Sj  vs  H2:  N  +  S2: 
GENERAL  THEORY  WITH  EQUAL  COST  ASSIGNMENTS 


Here  the  ternary  signal  detection  problem  is  considered  when  Ho:  N,  "noise  alone," 
is  an  alternative  state.  The  analysis  is  general  at  this  stage,  except  for  a  usually  noncritical 
assignment  of  decision  costs:  incorrect  decisions  are  penalized  equal  costs  C0  (>  0),  while 
correct  decisions  are  charged  zero  costs  (see  section  23.1(c)  of  reference  2).  However, 
there  may  be  unequal  a  priori  probabilities  of  signal  presence  (or  absence)  in  the  received 
data,  e.g.,  P2  =  (q,pi,P2)  ^  0,  q+pi+p2  =  1  for  the  present  case  of  two  possible  signal 
classes  (Si)  and  (S2).  (For  an  account  of  the  general  M-ary  theory,  M  >  2,  see  chapter  23 
of  reference  2.  See  also  references  14  and  15,  where  the  M-ary  detection  formulation  is 
used  for  signal  classification  in  Gaussian  noise.) 


2.1  THE  DECISION  PROCESS:  (M  *  2) 


With  the  cost  assignments  (see  pp.  1031  et  seq.  in  reference  2), 


C<0) 

c<k) 


C0;Z  =  1,2 
:C0;  /  =  k 


"failure" 


C(kk)  =0; /  =  k  =  l,2;  C<o)=0J 


and  from 
equation 
(23.7) 


Xio)  _  c<°)  —  C^0) ;  /  =  1,2  (=  M  here)' 
X(kk,sC(fck)-C^k)(<0);k  =  Z  =  l,2  ■; 
X(')  =  C(k,)-Ci,);/#k(=l,2) 


(2-1) 


thus, 

3l”=Cf-C®  =  0  1  A?«C®-CS"=-C.l 

x^cf’-c^c.J  a»=c®-cj*=-c.J  J 


with  the  thresholds 


(I)  _  A2  Ar,  A>2  #  J  A.2  A.} 


KL  = 

=  1 


;  W-C(0!)  ; 


=  1 . 


(2-3) 
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The  associated  decision  rules  (see  equations  (23.20a)  through  (23.20c)  of  reference  2))  are: 


decide  y0  if 

AU(V)S1;  (AkO); 

decide  y,  if 

A,(V)SA2(V)  and  A,(V)2>1; 

decide  y2  if 

A2(V)£At(V)  and  A2(V)£1, 

which  are  sketched  in  figure  2-1.  (The  simple  geometry  is  a  consequence  of  the  choice  of 
simple  costs  (see  equation  2-2).)  Figure  2- lb  shows  the  corresponding  decision  regions 
for 


xk  =log  Ak,  k  =  i,2  , 


(2-5) 


viz., 


decide  ya:  xt,x2  £  0  ;  decide  y,:  x,  £  x2 ;  Xj  <  0  ; 
decide  y2:  x2  £  x,  ;  x2  £  0  (see  equation  (7-4)) 


(2-6) 


(a) 

Figure  2-1.  Decision  Regions :  yb>7i,72>  and  Boundaries 
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The  equalities  in  the  decision  rule  relations  (see  equations  24  and  2-6)  define  the 
boundaries  of  the  decision  regions  Cyp,  j  =0,1,2  here2. 


2.2  BAYES  RISK 


As  seen  in  equation  (23.37)  of  reference  2,  specialized  to  M  =  2  (k  =  0:  noise 
only),  the  Bayes  Risk ,  or  Minimum  Average  Cost,  of  decisions  is  given  by 

RM=R«+qiw+  i>xH(  nr). . 


k-1 


kmi,lm0 


k»0 


(2-7a) 


where 


j-0 


(n  n=>-rv<fn 


(2-7b) 


for  the  present  ternary  case.  Using  the  cost  assignments  of  equation  (2-1)  yields  the 
ternary  Bayes  risk  here 


Rj  =qc0(a!“,'+o5’,‘)+c1,Jp1({p*1|")i+{pjl>‘)i)+p,((pi,!)')j  +{Pf’’)1)  • 


(2-8) 


with  9lo2  =  (p,  +  p2)C0 ,  q  +  Pi  +  p2  =1 .  For  a  check  consider  the  binary  case  M  =  2  , 
Si  vs  S2 ;  p,  +  p2  =  1  : 


•••R;=3t.!-c.[p,(i-{pr)1)+P:(i-{p?i')!)] 

=c0[(p,(p',>-)i)+(p2(p->-)J]. 


(2-9a) 


2  If  different  costs  are  assigned  to  the  various  "successful”  and  "unsuccessful"  (or  incorrect)  deciskn ,  then 
the  simple  decision  rules,  as  in  equation  24  above,  are  generalized  to  sums  of  likelihood  ratios,  rather  than 
ratios,  with  the  resulting  great  increase  in  the  analytic  complexity  of  the  decision  process  and  performance 
evaluation. 
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which  is  correct.  Also,  for  the  binary  "on-off1  cases  pj  =  q  ,  p2  =  p  ,  91  oM  =  pC0  ,so  that 


n-off 


(2.%) 


which  is  likewise  correct 


Rather  than  express  performance  in  terms  of  the  error  probabilities,  it  can  be  more 
convenient  to  use  the  (unconditional)  probabilities,  (P^) ,  of  correct  decisions,  when  it 

comes  to  their  evaluation.  With 


e-p,(i-(r), -(??•),).  =  PiPd ‘ ;  ' 
P<o’,'-Pi(l-(P<2)')!-(Pr):).  =  PiPo’-  ;  ■  . 

pr=q(i-«r-«r).  =q?r 


(2-10) 


where  P^*  are  the  conditional  probabilities  of  correct  decisions,  the  Bayes  risk  (see 
equation  (2-8))  can  be  written  more  compactly: 


(  2  > 

f  2  > 

Rm=2=C0 

i-ze 

o 

u 

II 

1 

j>0  > 

(2-11) 


Note  that  in  the  limiting  cases,  Sj,S2  -4  °°  ,  P(^*  -4 1 ,  and  therefore  P^*  -»  Pj ,  , 

2 

^pj  =  1  ;  therefore  R„=2  =0  ,  as  expected.  Conversely,  when  S,,S2  -4  0  ,  P^‘  =  Pj/2  ; 

0 

therefore  P^*  -4  pj/2  ,  and  therefore  R^=2  -4  C0/2  at  the  other  extreme;  for  these  cost 
assignments,  see  equation  (2-1). 
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2 . 3  CANONICAL  THRESHOLD  DETECTION  ALGORITHMS,  S,  AND  S2 

From  reference  16,  for  each  signal  class  Sj  andS2.  the  following  canonical3 
coherent  and  incoherent  threshold  detection  algorithms  exist  for  independent  sampling  of 
the  noise  field f  but  with  generalized  signals  and  noise. 


2.3.1  Coherent  Detection,  <s>  *0:  (k  =  1,2) 


logAgL=g^=A<°*'  +  Z,(»), 


coh 


(2- 12a) 


where 


A^slog  \i,-a 


■<k)J 

coh 


•  z> 

j 


(2- 12b) 


Here,  as  before,  k  =  1,2  refers  to  the  generalized  (and  normalized)  signals  while 
j  =  (m,n)  represents  a  space-time  sample.  A^>o, 


/(xJ)=^I°gw,(x|H.)|i_ii  =|i!w,(x|H„)dx-  |i,=^  ,  (2- 12c) 


with 


8,  ;  a«-)  =  )  =  1  ;  >P  =  I  +  4>G  (2-12d) 


3  "Canonical”  here  means  that  t be  forms  of  these  algorithms,  and  their  associated  performance  parameters, 
remain  invariant  of  the  particular  physical  application,  although,  of  course,  the  numerical  values  will 
depend  on  the  specific  problem  at  hand. 

4  This  latter  is  equivalent  to  treating  the  noise  field  as  spectrally  "white"  in  both  time  and  space.  As  noted 
in  section  1,  this  sampling  constraint  can  lead  to  serious  degradation  of  performance  vis-i-vis  algorithms 
which  property  include  at  least  the  dominant  first-order  correlations  of  the  noise  when  this  is  physically 
possible  (see  references  9  and  10).  This  situation  will  be  considered  in  a  subsequent  study. 
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for  the  normalized  signal,  with  actual  amplitude,  A*”’ :  m  =  sensor  or  "spatial 

sample"  in  the  received  field,  while  n  refers  to  ^ ,  the  n*  time  sample.  Here  I  =  intensity 

of  the  non-Gaussian  noise  component,  with  V'G  the  intensity  of  the  (frequently)  much 

smaller  background  Gaussian  noise  contribution.  Local  stationary  and  homogeneity  is 

assumed  over  the  bounded  sampling  region  ej  =  MN,  m  =  1.....M  ;  n  =  1.....N  .  For 

coherent  detection,  signal  epoch  is  known  and  equals  £<>  so  that  (s<0m))  *  0  ,  and  hence 

»  * 


2.3.2  Incoherent  Detection,  {s)=0 :  (k-  1,2) 
Similarly,  from  reference  16, 


(2- 13a) 


where 


Ar» = log  n,  -i<’ ;  zi(x)ik  *  XW) ; 


/;=— / 
dx 


;L(4’. {(/'+/>)’)  =|W  Wj(xjH„)  dx 
*,  '  /H»  VW1 


Also, 


fa2  *  ith  -  (m.ro')l  _ 

Vaom|a-n'|  Pjn-n'1  )k  ’  W1U1  "Vo'j  |k  ~ 

\aon  aoo'  / 

a 2 

v  °  y 

,  k  =  1,2  , 


p[nmn” }  ^  ,  with  (s(0m))  =  0  for  incoherent  detection 


The  narrowband  signals  considered  here  may  be  written 
s(am)  =  V2cos  [(o)o+codXta-e  +  ATm)-0n],  tn=nAt, 


(2- 13b) 


(2- 13c) 


(2- 14a) 
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in  which 


f0  =  "carrier"  or  central  frequency  of  the  signal  spectrum; 
o )d  (2jtfd)  =  Doppler,  often  consisting  of  a  determininistic  (f^) 

and  a  random  component  (f^); 
e  =  signal  epoch  =  e0 ,  known  for  coherent  observation  ; 
(usually)  uniformly  distributed  over  a  carrier  cycle, 
l/f0 ,  for  incoherent  observation; 

$■>  =  <j)(tn)  =  a  possible  phase  modulation:  for  example, 

PSK ,  where  4>B  =  rcn ,  t,  <  t  <  i,+I  ,  etc. 


where  one  or  more  of  the  signal  parameters  are  different,  each  k  =  1,2  .  The  quantity  Ax(mk) 
is  the  net  path  delay  from  the  m*  sensor  to  some  selected  reference  point  associated  with 
the  array,  viz.: 


(2- 14c) 


^  A 

where  iok  and  ioR  are  respectively  the  unit  vectors  in  the  direction  of  the  (incoming  kta) 
signal  wavefront  and  the  main  axis  of  the  beam  formed  by  summing  the  outputs  of  the 
(receiving)  sensors  (lm) ;  rm  is  the  vector  distance  of  the  m*  element  of  the  array  to  the 
reference  point;  Co  and  are  respectively  the  speed  of  (group)  propagation  and  the 
(central)  wavelength  of  the  (kih)  signal  in  the  medium. 


With  incoherent  detection  [(e  *  e0) ,  (s)t  =  0  ,  or  (0)  =  0]  ,  the  signal's 
autocovariance  (9j0j<)  ,  instead  of  the  waveforms,  are  dealt  with  (see  equation  (2- 13c)). 


The  carrier  or  rapidly  varying  portion  of  (OjOj-)  becomes 

cosjo)„(n  -  n')At  +  k,(i.  -  ?„)•  Ar„.] 


Pf*?,  =  1  exp: 


Ao)d(n-n')At] 

2 


\ 


k  = 

v  ° 


,(2-15) 


where  Acod  is  a  measure  of  the  "Doppler  spread,"  if  any,  assumed  to  be  Gaussian,  while 
the  [<t>n]  are  independent  over  n  here,  and  Arm  m,  =  rm  -  rm, .  Clearly,  directing  the  beam 
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so  that  ioR  ->  iok  maximizes  p,  as  expected,  while  uncertainty  as  to  Doppler  (Ao)d  >  0)  and 
signal  direction  (ic)  degrades  p  and  hence  performance  (see  references  12  and  16). 

Since  appropriate  delays  and  weightings  for  each  sensor  output  can  be  inserted 
here,  there  is  an  opportunity  for  adaptive  beam-forming.  The  space-time  algorithms  in 
equations  (2-12)  and  (2-13)  provide  (threshold)  optimal  "rules"  for  so  doing,  to  the  extent 
permitted  by  learning  during  data  acquisition  (s  "adaptive  processing")  and/or  a  priori 
knowledge  of  the  signal  and  noise  parameters  and  statistics,  e.g.,  a* ,  a2 ,  e0 ,  w,(x|H0) , 
etc.,  required  to  implement  these  algorithms.  This  adaptive  beam-forming  and  the 
associated  temporal  processing  (tn)  is  called  decision  matched-field  processing  (see 
footnote  1),  which  is  optimal  from  the  point  of  view  of  the  general  criterion  of  decision 
theory  and  Bayes  risk.  This  is  more  general  than  the  usual  criterion  of  matched-field 
processing  for  estimation,  or  detection,  based  on  maximizing  array  output  on  the  basis  of 
signal-to-noise  ratio5. 

In  many  instances,  the  receiving  beam  is  preset  or  preformed  with  no  option  of 
adaptive  beam  forming.  Then,  formally,  the  spatial  index  in  the  above  algorithms  is 
dropped:  j  — »  n  ,  J  -4  N  ,  and  M  =  1 ,  with  Ax^  -4  0  ;  only  time  processing  is  permitted. 
This  is  the  usual  situation  in  most  electromagnetic  (EM)  telecommunication  applications, 
where  array  L  <  wavelength,  but  not  in  underwater  acoustics,  where  L  »  \0 ,  and  thus 
the  structure  of  the  noise  and  signal  fields  permits  its  exploitation  for  improved 
performance. 


2.4  PERFORMANCE  PROBABILITIES:  TERNARY  (M  =  2) 

At  this  point  there  is  a  transform  to  the  new  random  variables: 

Z=U-A ; A  =  (Ak) ,  U  =  (log  A(k))  ;  Z  =  (Zk)  ,  k  =  1,2  .  (2-16) 


5  When  dealing  with  parameter  estimation,  the  appropriate  decision-theoretic  extensions  must  be  used  (see 
reference  2,  chapter  21,  for  example)  based  on  maximum  livelihood  estimation  (MLE),  or  least  mean-square 
error  (LMSE)  criteria,  which  follow,  in  turn,  from  the  Bayes  risk  theory  by  appropriate  choice  of  cost 
function  (see  section  DC  and  the  appendix  of  reference  16,  for  example).  (See  also  reference  17  for  an 
alternative  approach  in  underwater  acoustics  far  estimating  source  location.) 
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Next,  use  equations  (5)  and  (6)  of  problem  23.2  in  reference  2,  to  write  the  various 
conditional  error  probabilities  (M  =  2)6 : 


(2- 17a) 
(2- 17b) 
(2- 17c) 


These  relations  are  also  readily  established  from  figure  2.  lb  from  the  definitions  of  the 
error  probabilities.  Here  Q2  and  P2k)  are  the  PDFs  of  Z,  (see  equation  (2-16)),  to  be 
determined  presently  for  the  canonical  threshold  detection  algorithms  in  equations  (2-12) 
and  (2-13)  and  in  section  3. 

The  forms  of  the  probabilities  of  correct  detection  follow  in  similar  fashion  from 
equation  (2-10): 


(2- 18a) 
(2- 18b) 

and 

(2- 19a) 
(2- 19b) 


’£'■  =  P  JdU,  [  JdU, + )  du, W(U) = p,  j  dU,  ]dU2P*’»(U) , 

o  v-  o  /  o 


Zt+Ai  *Aj 


=  p,Pg»-  =  p,/dZ,  JdZ^fZ). 


Pd**  =  qJdU,  JdU,Q,(U) ,  U  =  (Ut,U2) (see equation  (2-16)) , 

-A,  -A, 

=  qpr=qJ<iZ,  JdZ,Q,(Z), 


6  Note  the  (inadvertent)  omissions  of  -Aj  and  -A2  in  the  limits  of  <£* y*>  and  <P(21)*>  in  equation  (A.l-14c) 
of  reference  3  and  correspondingly,  in  p(Jj*  and  p^*  (see  equations  (A.l-16b)  and  (A.l-17b)  of  reference  3). 
The  consequences  of  this  are  noted  in  appendix  B  here. 


2-9 


Similarly, 


•• 

f  0  u,  \ 

—  Uj 

P™'  =  pJdU, 

JdU,+  JdUi 

P<2,,(U)  =  pJ<lUI  JdU,Pf(U), 

0 

\-~  o  J 

«•  Zj+Aj'Aj 

0  — 

=  pA2,‘  = 

Pj|dZ,  JdZ,P™(Z). 

_ '±i _ z _ 

(2-20a) 

(2-20b) 


As  measures  of  performance,  the  P^*  or  p^*  (where  j  =  0,1,2)  can  be  used  directly  or  in 
conjunction  with  equation  (2-11)  for  the  Bayes  risk. 
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3.  EVALUATION  OF  THE  PERFORMANCE  PROBABILITIES: 

TERNARY  (M  =  2) 


The  first  task  is  to  obtain  the  PDFs  Q2  and  P2k) ,  of  the  (here)  threshold  detection 
algorithms  in  equations  (2-12)  and  (2-13),  and  then  to  determine  the  various  decision 
probabilities  in  equations  (2-17)  through  (2-20)  associated  with  this  canonically  formulated 
ternary  detection  problem.  Two  critical  features  permit  comparatively  simple  analytic 
results,  expressed  at  worst  by  convergent  series  of  known,  tabulated  functions.  These 
features  are:  (1)  the  asymptotic  normality  (A-N)  of  the  algorithms  g^  in  equations 
(2-12)  and  (2-13)  by  the  Central  Limit  Theorem;  (2)  this  ternary  situation,  which  involves 
only  two  (M  =  2)  signal  classes  in  addition  to  noise  alone,  which  fact  is  reflected  in  the 
double,  rather  then  triple,  integrals  representing  the  desired  decision  probabilities. 


3.1  THE  NORMALIZED  PDFs  OF  Z 


The  A-N  nature  of  g* ,  because  of  the  weak-signal  constraint,  allows  the  needed 
PDF  of  the  Zs,  now  in  normalized  form,  to  be  written  as 


rxn,  (Z' -2zAP  +  z’) 


,(T1  =  0,1,2) 


qJ(Z),T1  =  0, 
_P;k)(z) ,  fl  =  k  =  l,2  ,J 


=  P2(Z). 


(3-la) 


See  section  7.2  in  reference  2,  where  specifically 

_zk-(zk), 


k  =  1,2 

(n =o,i,2) 


:  i 


»  <Zk)n*(Zk(x))Hn  ;  (  >„  s  JWl(x|H„)dx ; 

wl  — 

p  =  p?  2  =  [(Z,  -  <Z, ),  XZ2  -  <Z,  >2  ; 

o™'.var,,<Z1=[(zl-{Zt),)J] 


(3.1b) 
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Note  that  Pj\z)  =  p{22)(z)  in  this  normalized form  (see  equation  (3-la)). 

In  sections  3.1.1  and  3.1.2,  the  various  moments  that  are  summarized  in  equation 
(3-lb)  are  evaluated  using  the  results  of  appendixes  A-l  and  A-2  of  reference  18. 


3.1.1  Coherent  Detection 


<Z«>. 

<L 


=  0  (see  equation  (A.1-  13b)) ; 

=  L<2)]£{e5k))  ■  (see  equation  (A.2  - 14)) ; 

j 


H, 


0  * 


(°i02)«*Pi°2  ”  (ZjZj)jjo  =L(2>^(0(jI)^0f))-(tT,Cf2)cohpcoh 


since  (Vr)H#  ■  . 


(3-2) 


Similarly,  (t|  =  1,2) : 

(Zk)w  =  L(2)X(0jk))2  +  Op)  (**  equation  (A-2  -  3) ,  =  o2k.rob  ; 
j 

c(kn)J  =  L(2)X{0jIC))2  +  o(0<)  (see  equation  (A.2  - 3) ,  = 
and  0,C2pit  =  LU)X(ejI>X0J2>)+ (K0*) 4(a»°aLP«*»  • 


3.1.2  Incoherent  Detection 


(Zk)0  =  0  (see  equation  (A.2  -  25) ; 
CTk°iac  =  equation  (2  -  13b)  exactly ; 


=  (°iCTzLPfcc .  exactly. 


(3-4) 
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For  the  alternative  hypotheses  (k=l,2) , 

(Zk)]  2  =  =  equation  (2.13b) ;  =  CT2k_Uc  (see  pp.  226  -  232) ; 

(oS’,)o(2n)p|k2))jac  =  (a^J^p*.  (see  equation  (3  -  4)) 

From  equations  (2-12d)  and  (2-13b), 

A^elogp,-^^  ;  A^slogp,--^  ;k  =  l,2,  (3-6) 

using  equations  (3-2)  through  (3-5)  for  the  detailed  structures. 

Accordingly,  to  summarize,  equation  (3-16)  can  be  compactly  expressed  as 

(3-7a) 


and 


(3-7b) 


In  the  subsequent  analysis,  the  subscripts  "coh"  and  "inc"  on  zk ,  crk  ,  etc.  will  be  dropped 
unless  this  results  in  ambiguity. 
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3.2  THE  PERFORMANCE  PROBABILITIES:  (M  =  2) 


Using  the  transformations  of  equations  (3- lb)  on  the  unnormalized  (Gaussian) 
PDFs  for  Zi  and  Z2  (see  equation  (7.9)  of  reference  2),  gives  the  equation  (3- la),  as  noted 
above,  and  allows  Q2(Z)  and  P(2k)(Z)  in  equations  (2-7)  through  (2-20)  to  be  replaced  by 
q2(z)  and  p2k)(z) ,  respectively  {see  equation  3-1)).  The  results  for  equations  (2-17). 
through  (2-20)  now  take  the  explicit  forms: 


(3-8a) 

(3-8b) 

(3-8c) 


Note  that  since  p2k)(z)  =  p2(z)  in  equation  (3-1),  the  subscripts  1  and  2  can  be  interchanged 
in  the  limits  of  all  the  error  probability  integrals  to  obtain  the  (p{  }*)s  in  equation  (3-8). 

The  same  can  be  done  for  equations  (2-18)  through  (2-20)  in  normalized  form 
(without  the  respective  a  priori  probability  factors): 


(3-9a) 

(3-9b) 


-A,/«,  —  Aj/Oj  ..  (oi*i+A,-A,+erf-ol)/<J2 

P(D0)‘=  Jdz.  |dz2q(22)(*);  Pd*  Jdzi  |dz2p2)(z) ; 

-  (aiij+Aj-A^ol-ofJ/o, 

Pd}*  =  Jdz2  J dz,p(2J)(z) ;  and  p2k)(z)  =  p2(z)  (see  equation  (3  - 1)) . 

_ (-Ai-gn)/g; _ 


m  (0|Z,  ■fA1-Aj)/Oj  (OiZj+A2-A,)/Oi 

«!0>‘ePo),=  Jdzi  }dz2q2(z);  a(20)‘  =P(20)*  =  Jdz2  Jdz,q2(z); 

-A,/<»|  —  -Aj/Oj  — 

(-A.-O?)/®,  (-Aj-aj  )/oj  „  (fljZj + Aj  -  A,  +«!  -<f  *  )/<J, 

(Po)*)=  Jdz!  Jdz2p(21)(z);  ($Y)  =  Jdz2  Jdz,p(2,)(z) ; 


(-A,-o?)/a,  (-a2-o’)/<j2 


-Aj/Oj 


(Po2)‘}=  Jdzi  J dz2p22)(z) ;  (p{2)‘)=  Jdz, 

_ —  _  —  _  (-A,)/o, 


(o.z.+A.-Aj+of-olJ/o, 


Jdz2p(22)(z). 


3.2. 1  Probabilities  of  Correct  Detection 
Next, 

q2(*) .  p2k)(z)  =w2(z]Hn)  equation  (3-1) , 

=  w2(x,,x2)  (see  equation  (A-l)  in  reference  19) 


(3-10) 
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and  all  have  the  same  analytic  form  (see  equation  (A-l)),  so  the  basic  result  of  I  =  equation 
(A-6c),  with  equation  (A-8b)  and  special  cases,  is  to  evaluate  all  the  decision  probabilities 
in  equations  (3-8)  and  (3-9).  These  results  are  specifically: 


Explicitly,  the  following  series  may  be  used  (for  |p|  <  I/V2)7 


(3-1  lb) 


where  Ic(a;v)  s  (l/r(v))J yv_1  e'y  dy ,  v  >  0  :  incomplete  F-  function  (see  reference  22) 

0 

and  <J>(m)(x)  =  dm(e"xl/2/V27t)  dxm  is  tabulated  also  (see  reference  19).  Here  0(x)  = 

X 

(2/  Vtt )  J  e'*1  dt  =  erf  x  is  the  familiar  error  function  of  references  20  and  21. 

0 


7  For  |pf  <  1/V2  ,  another  series  is  needed.  This  is  developed  in  reference  13,  allowing  a  complete 
numerical  evaluation  for  -1  S  p  £  1  . 


Further  inspection  of  equation  (3-1)  in  equation  (3-9a)  shows  that  if  p  =  0  , 
then  w2  =  w,  (zt  )w,  (z2 ) ,  i.e. ,  w2  is  symmetrical,  and  so  equation  (3- 1 1  b)  reduces  to  the 
much  simpler  relation 


pa.4{H^)H 


-r2- 

V2o2 


(3-12) 


whereas,  if  p  *  0  ,  then  p(D0)‘  is  unsymmetrical  (in  Aj/Oj  and  A2/o2),  as  expected,  since  w2 
is  also  unsymmetrical.  In  the  limits  p-»±l ,  w2(z)=  Wj(z,)8(z2  Tz1)  =  w,(z2)8(z1  +  z2),  and 
equation  (3-9a)  must  be  evaluated  directly  to  get 


p»’‘L ;  B-  ■  mi"(^  •  ■ 

(where  min(-5,-2)=-5  ,  etc.).  (For  the  cases  ,  <J2  — >  0  ,  see  equation  (3-17).) 

The  conditional  probabilities  p(Dk)*  (in  equation  (3.9b))  are  found  in  the  same  way. 
Equation  (A-6c)  becomes 


(3-13) 


Pif=- 


2V2tt 


-P 


A,  -  A2  +  of  -  a2 


V2(1-p2)  *  cr2V2(1~p2) 


-J 


-A.-q? 


-P 


A, -A2-Hg?-q2 


V2(l“P2)  ’  02^2(!-  p2) 


+J 


-p 

—  J 

1 

1 

> 

1 

Q 

— 

P 

J 

V2(i-p")  ' 

(3- 14a) 
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which  reduces  (for  |pj  <  1/V2 )  (see  footnote  7)  to 


(3- 14b) 


Note  that  for  (<J2  <  «) ,  (A,  +  <jf  )/o,  -» a,/ 2  — >  «> ,  etc.,  Ic  ->  1  ,  and  p^*  ->  1,  as 

required:  Si  is  surely  detected  (conditionally  on  being  present).  In  this  limit  (o,  — » «») , 

«• 

equation  (3-9b)  becomes  also  p^*  =  J  J  p^(z)  dz  =  1 ,  as  a  further  check.  Since  p(Dk)*  is 

symmetrical,  e.g.,  1  ->  2  ,  2  ->  1  in  equation  (3-9b),  it  can  be  written  explicitly  (|p|  <  l/  V2)  (see 
footnote  7): 


Again,  as  expected,  p^*  1  as  c2  -»  «  . 


3.2.1. 1  Special  Cases.  When  p  =  0  ,  then  p(Dk)\  k  =  1,2  ,  does  not  reduce  to  the 
symmetrical  form  of  equation  (3-12)  for  pj)0)*,  because  of  the  asymmetries  in  signal  level; 
o,  *  o2  .  Furthermore,  note  that  for  O) ,  o2  -»  0  ,  since  4>(m)(±<»)  =  0  ,  m  >0  : 
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(3- 16a) 


lim  J_ 

Pd  a->0  4 


■*K® 


1  +  0 


log 


t 

Pi 


V2o" 


Several  cases  can  be  distinguished: 

Pi, 2 


Hu  >  1 


>  1 : 
P2.1 

^<1: 

P2.1 

M-i  2  <  1 ;  —  <=>  1 

P2.1 

Hi  -  h2  :  Pi  =  p2 ;  Pi  ,2  >  q : 


P(Dk)* 


P1.2  <  q 


1 


Pi  =  p2  =  q  =  3 ;  h,  =  h2  =  1 : 


,  j,k  =  1,2  ;  j*  k 


=  1; 
=  0; 

=  0; 

_J_ 

"2  1 

=  0; 

1 


;  k  =  l,2 


For  the  noise-alone  case  in  equation  (3- 1  lb),  since  <J>(m)(±°o)  =  0  ,  m  >  0 
lim  1 


pS'  = 


ou-+0  4 


1  +  0 


^-logHi^ 


Via 


1  n 


1+0 


^-logHz^ 


Via 


2 


(3- 16b) 


(3- 17a) 
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and  therefore 


M-1,2  >  ^  * 

P(D0)‘  =0; 

Pu  < 1 : 

=  1 ; 

Pj.2  >  ^  *  P2.I  <  ^  • 

=  0; 

p,  =  p2  >  1 : 

=  0; 

p,  =  p2  <  1  : 

=  1; 

,  1 

1 

Pi  *  Pa  *  1  .  ”  P1.2  =9  =  3  : 

=  4 

Pl,2  ^  »  P-1,2  =  ^  * 

_  1 
”  2 

(3- 17b) 


3.2.2  Probabilities  of  Error 

The  error  probabilities  of  equation  (3-8)  can  be  calculated  in  similar  fashion  with  the 
help  of  equation  (A-6c).  The  general  results  for  equation  (3-8a)  are 


««»■  - ,  1 
1  2V2S 


--P 


^  Aj  A2 

V2(i-p2)  ’  <W2( >-p 2) 


-J 


-A, 


Oj 

— — p 

o. 


_  A, -A, 

V2(i-p2)  ’  o2^2(l-p2) 


+J 


^M2)  ’ 


"A, 

-P 

a/2(‘-p!) 

Upj0)*.  (3- 18a) 
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The  series  form  (for  <  l/V2)  of  equation  (3- 18a)  is 


a 


(0)*  _ 


i-el 


V2 a. 


2  ^-1  1 


- 

/ 

1  +  0 

i 

TJ. 

O-j 
V  2 

V2 

p)l 

J 

^1-P: 

rf -2^1' \b(,-,>f A; 


'y—r  rai-p- 


1-1, 


A:  n  +  1 


.2°i  2  . 

|p|  <  l/V2  (see  footnote  7) . 


(3- 18b 


Note  that 


a: 


(0)- 


0(-»-  =  0  .  while  a(20)* 

o,<- 

-I 

o,<-  ^ 

(  _2  V 

-lOgm+^L 

l-e 

V2CT, 

. 

l  A 

(3- 18c 


For  a^0)’(s  #0)*) ,  the  subscripts  1  -*  2  and  2->  1  are  exchanged  in  equations  (3-1 8a) 
and  (3-18b). 


Next,  for  equation  (3-8b), 


-A,  -  a? 

-P 

-Ji 


-P 


A2  -  a] 


W?)  * 


UpS*.  (3*  19a' 


+J 


-A.-g' 

-  P 

-  J 

P 

<*1 

#-pV  J 

J 

V2(*-ps) 

3- 


I 
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|p|  <  l/V2  (see  footnote  7) 


Since  p(2y)  =  p\2)  =  p2(z)  (see  equation  (3-1)),  here  (f}(02)*)  =  (Z^*)  a&o  (/rom  equations 
(3-8b)  and  (3-8c))*. 

Similarly,  for  and  (p|2)*) 


8  In  fact,  the  equality  '  J  *  '  j  here  (see  equation  (3-8b),  with  P2  (7)  =  P2(z)  •  each  lc, 

allows  the  interesting  equivalence: 

•  “t)-' ■ 


1 


and 


For  (/>j2)*)  ,  the  subscripts  l—>  2  and  2-4 1  in  equation  (3-20b)  can  be  interchanged 
(compare  equations  (3 -8b)  vs  (3 -8c).  Observe  that  as  Oj  or  o2  -» «> ,  these  (conditional)  error 
probabilities  (if  signal  1  or  2  is  actually  present)  all  vanish  as  expected  from  physical 
considerations  (see  section  4). 


3.2.2.1  Special  Cases.  In  addition,  as  O,  or  a2  -4  0  ,  (see  equations  (3-16)  and  (3-17)), 
there  are  a  variety  of  special  results,  depending  on  the  a  priori  probabilities  q  ,  p_  ,  and  p2 . 
These  are  determined  from  equations  (3-18)  through  (3-20),  which  become,  when 
(<y  =  o1,a2)-40: 


(3 -2  la) 

<Tj  ,  Cj  =  0  -4  0  . 

(3 -2  lb) 
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As  before  (see  equation  (3-16b),  several  cases  can  be  distinguished: 


Pi. 2  >  i ; 

£!i>1:  0W 

p2.l 

=  1; 

Pu<l; 

P2.1 

=  0; 

Pu  <  i ; 

^2.<=>1: 

P2.1 

11  11 

to  1  —  O 

>  c-»0 . 

(3-2  lc) 

Pi=p2  pi  =  p2 ; 

Pu>q: 

Pu<q- 

=  1; 

l 

Pi  ~  p2  -  q  ~  ^  ’  *** 

Pi=P2=l: 

4 

(Note  that,  in  fact,  equations  (3-21a)  and  (3-21b)  are  the  same  as  equations  (3-16a)  and 
(3-16b).) 


In  a  similar  way,  when  o, ,  <J2  0  ,  then  equation  (3-19b)  reduces  to 


(C) 


lim  1 
a-* 0  4 


1  + 


(3-22a) 


The  various  cases  here  are  therefore 


Pu  > 1  • 

(fTHfl r)  =0; 

Pu<1: 

»1; 

|ii  >1  •  M-<1 : 

=  0; 

m>i,n2*ii 

=  0 

t 

M-2  >  1 ,  M-,  =  0J 

=  0 

P,  =  >  1] 

=o: 

Pi=P2<U 

=  1 

-  1 

Pi =  p2  =  q » Pi =  p2  = 

1 :  =- 

4 

Pi = 1 ;  p2  <  1 ,  Pi  =  q  ] 

L.  =11 

P2  =  1 ;  Pi  <  1 »  p2  =  qj 

2J 

(3-22b) 
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Finally,  for  and  (p|2)*)  ,  from  equation  (3-20b),  (Oj  2  -»  0) 

l0ffe)l 


(or) 


lim  i 

=>  „  — < 
o  ->  0  4 


V2o 


Al. 


and 


lim 

o-»0 


(3-23a) 


(3-23b) 


Again,  note  that  since  equations  (3-23a)  and  (3-23b)  ate  the  same  as  equations  (3-21b)  and 
(3-21a),  respectively,  the  various  limiting  cases  of  equation  (3-21c)  apply  directly  here  as 
well. 
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4.  ADDITIONAL  SPECIAL  CASES 


Further  special  cases  of  interest  arise  when  one  of  the  (processed)  signals  is 
indefinitely  large,  e.g.,  when  c,  and  c2  <  «> ,  or  o2  < «» .  See  equations  (3-1  lb), 
(3- 14b),  and  (3-15)  for  the  p(D}*s  and  equations  (3-18b),  (3-19b)  and  (3-20b)  for  the 
following  results. 


4.1  (<Tt  2  ->  oo  ;  cr2  1  <  oo) 


For  the  probabilities  of  correct  detection. 


(4 -la) 
(4 -lb) 
(4  -  2a) 
(4 -2b) 


Because  of  the  signal  level  asymmetries  Oj  *  <x2 ,  c,/o2  -»  °° ,  and  p  *  ±1 ,  p^*  may  be 
expected  to  be  nonvanishing,  even  when  one  of  the  signals  is  effectively  infinite:  there  is 
always  a  nonzero,  non-unity  probability  of  falsely  deciding  signal  #1  or  #2,  as  is  evident 
from  equations  (4-la)  and  (4-lb).  On  the  other  hand,  as  in  equations  (4-2a)  and  (4-2b), 
PDk)  L  =  1  •  PdH  =0  ,  j  *  k ,  and  ok  < «»  ,  as  expected;  very  («>)  strong  signals 

KJj 

are  correctly  decided,  probability  1,  with  the  alternative  p(D  =  0  ,  by  exclusion. 
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In  a  similar  way,  for  the  error  probabilities: 


'Oj  <-  'Og  < 


.  /v<0)* 


a20)*|«f,-»-  =  0  ;  a^0)*L-»-  - 

'o,<-  ,02<~ 


r  ~  -a,  * 

> 

l-e-p- 

L 

• 

* 

r  -a,  " 

i-e-7^2- 

L  ^<*2. 

. 

and,  from  equation  (3- 19b), 

(r|,.-=»=r);(fL-=« 

'Oa<- 

and,  from  equation  (3-20b), 

(P«%,~  =  0;(rl~=0; 


and 


riaz<^»  'aj  <•• 


=o 


(4 -3a) 
(4 -3b) 


(4-4) 


(4-5a) 


(4-5b) 


Equations  (4-4)  and  (4-5)  combined  in  equation  (2-10)  yield  equations  (4-2a)  and  (4-2b), 
as  required.  Thus,  whenever  the  true  state  is  "noise  alone”  H0:  N ,  as  in  equations  (4-3a) 

and  (4-3b),  and  the  non-infinite  signal  is  decided,  the  associated  decision  error  probability 
is  nonvanishing.  Whenever  the  true  state  is  Ht:  S,  +  N  or  H2 :  S2  +  N ,  and  either  signal 
is  infinite,  the  error  probabilities  vanish.  Both  effects  are  expected:  in  the  former  case, 
there  is  always  a  nonzero  probability  of  incorrect  decisions  because  one  of  the  signals  is 
finite;  in  the  latter  case,  the  infinite  signal  ensures  a  correct  decision  or  "signal,"  and  hence 
zero  probability  of  the  alternative  decisions. 


4.2  THE  CASES  p  -  ±1 

Here,  as  in  equation  (3-12),  w2(z)  =  w,(z,)S(z2  Tz,)  =  w1(z2)8(zJ  T  z2) ,  so  that 
evaluating  equation  (3-9)  directly  gives  (the  earlier  result  of  equation  (3-13)): 
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rf-L  =il1+e(B-.)]  i B-  *  mi”(^  •  • 

4[,+(fe)L :  "Kj^L  ■ 


(4*6a) 


along  with 


D<*)*|  =  ©  fti ft?,,* £l _  ©f 

U±i  ^  V2(c2-o,)  J  \ 


-A  ,-qf 

V2o, 


2  A 


=  0, 


.  ~Ai  ~ai  At  -  A2  +  a,  -  q2 

>  'o 

CT j  (Jj  ""  0"| 

-Ai-Q?  >  At  -  A2  +q?  -q2 
°1  °2  ~  CT1 


(4 -6b) 


(4  -6c) 


and 


-Ja2  -A.+qj-g; 


— «i»— a— — 

V2(o, -0l) 


“9 


=  0. 


-A;-o]  eA,-A,+oj-o? 

^2  CTj  **“  0^2 

-A,-<^;  A,— A,+o»— or? 


CT,  O, 


.  (4-6d) 


In  a  similar  fashion,  from  equation  (3-8): 


-A, 


V2o, 


i  J 


a 


(0)- 


p»±l 


=  e 


=o, 

A2  ~  A, 

V2(o,-o2) 


.  At  A2  A,  # 

»  ^  t 

°2  “  ^1 

A|  ~  A;  ^  -A, 
°2  —  °1 


•€N 


f  _A  N 
a2 


V2a 


2  / 


=  0, 


A2  A,  ^  A2 
o, -a2  q2 

A2  ~  A;  ^  "A; 

o,  —  Oi  <T, 


Ai  a,  A2  < 


(4-7a) 


(4-7b) 


(4-7c) 
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;  (4-7(1) 


Ja2-a 

ef  ”A0 

V2(a,  -a2)  _ 

0, 


A2-At  +  CT*-o?  >  -A2 
(°i-°2)  o2 
a2-a ,+q^-g?  ,-a2 


(^-o2) 


and 


=  0, 


A,-A24-CTf-(^  ^  -A, 
(°2  ~  °1  ) 


/r<2)*M  _  n[ Ai  ~  A2  +  °i  - q2 1  J  "A1  V  A,  -  A2  +  of -c\  -A,  . 

V'  K«-*[  V^-c.)  J  (O,-*,)  "5“’ 


.  (4-7e) 


4.3  THE  CASE  p  =  0 

When  p  =  0  ,  the  situation  of  orthogonal  signals  exists  in  the  coherent  cases: 

S0  -  Vlaocos  a)0(tn  -ec) ;  8(„2)  =  V2a0sin  m0(tD  -e0) ;  ec  =  (tc/2g)0  .  3n/2co0) 
or  (0,tc/g>o)  for  (maximum)  amplitudes  for  proper  choice  of  epoch  e0  at  the  receiver,  e.g., 

SPS?  =  Pcoh  =  0  •  With  incoherent  detection  (0U  =  o) ,  (Q^e^V)  =  8nn. ,  etc.,  for 

p  =  0  .  In  any  case,  when  p  =  0  ,  the  general  results  in  equations  (3-1  lb),  (3-14b),  and 
(3-15)  are  directly  modified  on  setting  p  — >0  therein,  as  in  equation  (3-12).  Equations 
(3-14)  and  (3-15)  do  not  reduce  further,  unless  o,  =  o2  ,  in  which  case 


(4-8a) 


(4-8b) 


with  Pp here  given  by  equation  (3-12). 
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4.4  REDUCTION  TO  THE  BINARY  CASES:  Hu  vs  H2  AND  Hj  vs  H© 


The  more  general  results  (in  section  3.2)  readily  reduce  to  the  binary  cases  upon 
dropping  the  inappropriate  error  probabilities  and  regarding  q  as  q  -»  0  ,  e.g.,  log  (p,  ,  p2) 
— »  «>  with  p,  +  p2  =  1 ,  etc.,  now.  Thus,  for  binary  signals: 


H}:  S,  +  NvsH2:  S2  +  N: 


n(0*  _  1 
p°  -2 


f  -2  2 

gl  g2 


1  +  0 


+  log  p12 


;  log  1 1, 2  =  log 


Pl 

P2. 


and 


n(2)*=I 

Po  2 


1  +  0 


(£z°L+ 

2 


log  m , 


/J 


;  log  p2i=iogf— 
VPi> 


(4-9a) 


(4-9b) 


Equations  (4-9a)  and  (4-9b)  are  alternate  forms  to  the  more  familiar  results  (see  section 
20.4-5,  problem  20.12  of  reference  2),  which  latter  are  given  in  terms  of  log(A2/A,)  = 
log  A,  2  rather  than  log  A2  vs  logAj ,  as  is  done  here. 


Also  in  the  "on-off  cases  H0:  N  vs  Ht:  S,  +  N ,  now  p  =  0  ,  <r2  -+  0  ,  so  that 
equations  (3- 18b)  and  (3- 14a)  reduce  to 


1-0 


_ log  Hi 

2V2  V2o, 


J) 


(sa‘),  (K  =  1  here) ; 


(4- 10a) 


and 


(4-1  Ob) 


which  are  well-known  results  (K=l  here)  (see  equations  (6.6)  and  (6.7a)  in  reference  12, 
for  example). 


4-S/4-6 
Reverse  Blank 


5.  CONCLUSIONS 


The  principal  new  features  of  this  preliminary  investigation  are  the  explicit 
development  of  threshold  Ternary  Detection  Theory  for  the  two-signal  cases:  K0:  N  vs 
H,:  St  +  N  vs  H2:  S2  +  N  .when  the  noise  is  (additive)  non-Gaussian  and  when  the 
signals  are  themselves  unrestricted.  In  short,  this  is  a  canonical  theory,  invariant  in  the 
form  of  the  analytic  results.  Both  threshold  algorithms  and  detection  performance  are 
derived.  Special  cases  of  interest  are  then  readily  obtained  by  specializing  these  general 
results,  as  sections  3  and  4  attest 

However,  the  treatment  is  not  completely  general  in  that 

1 .  Threshold  signals,  not  signals  at  all  levels  of  intensity,  are  considered; 

2 .  A  restricted,  though  important  class  of  cost  assignments,  which  determines 
the  decision  regions,  is  employed  (see  section  2.1); 

3 .  Independent  noise  samples  are  assumed,  in  both  space  and  time; 

4 .  The  analysis  is  confined  to  the  two-signal  cases  (M  =  2),  with  "noise  alone"  as 
the  third  alternative  (see  section  2); 

5.  The  noise  is  ambient  and  thus  signal-independent;  and 

6.  For  numerical  results,  series  valid  only  for  |pj  <  1/ V2  are  obtained  here. 
Reference  13  contains  the  complete  range  for  -1  £  p  £  1 . 

Under  these  constraints  only  three  parameters  (other  than  the  a  priori  probabilities,  q,  pi, 
and  p2)  are  required  for  a  complete  analysis:  the  signal  detection  parameters  erf ,  erf ,  from 
equations  (3-2)  and  (3-3),  and  the  correlation  coefficient  p,  2(=  p)  relating  to  the  two 
detection  algorithms,  Z,  and  Z2 ,  from  equation  (2-16),  as  shown  by  the  results  of 
section  3. 
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Restriction  3  can  be  removed,  as  a  result  of  recent  work  (references  9  and  10).  The 
extension  to  the  full  ternary  case  H,:  S,  +  NvsH2:  S2  +  N  vs  H3:  S3  +  N  ,  item  4  under 
restriction  2  above,  appears  analytically  possible  though  tedious.  Restriction  5  can  also  be 
lifted,  as  noted  in  reference  1 1  and  references  therein.  On  the  other  hand,  analytical  results  for 
general  signal  levels,  and/or  general  cost  assignments  are  not  achievable,  and  computational 
procedures  must  then  be  invoked  for  specific  cases.  Accordingly,  next  steps  in  the 
development  of  the  ternary  detection  theory  here  are: 

1 .  Extension  of  (M  =  2)  threshold  cases  to  include  noise  correlation  (see 
references  9  and  10); 

2.  Numerical  evaluation  of  the  general  results  for  (M  =  2)  obtained  in  section  3, 
for  both  the  "correct”  detection  and  error  probabilities,  for  a  usefully  broad 
spectrum  of  values  for  the  three  parameters  involved,  viz.,  erf  ,  o2  •  311(1  P 
(see  section  3.1).  This  is  done  in  reference  13; 

3 .  Extension  of  die  present  results  to  signal-dependent  noise,  which  is 
particularly  important  in  active  regimes,  where  reverberation  is  dominant;  and 

4.  Formal  extension  to  the  full  threshold  signal  case  (M  =  3),  at  least  to  the  point 
of  identifying  ways  of  analyuc  evaluation  of  the  triple  integrals  (with  3rd-order 
Gaussian  integrands)  specifying  the  performance  probabilities. 

Finally,  there  remains  the  application  of  these  general  results  to  specific  problems  in 
sonar  and  related  "active"  as  well  as  "passive"  situations. 
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APPENDIX  A 

EVALUATION  OF  TERNARY  PROBABILITY  INTEGRALS 


In  the  analysis  developed  in  the  preceding  text  for  the  important  ternary  case  H0:  N 
vs  fy:  Sj  +  N  vs  H2:  S2  +N  ,the  general  probability  integral  is  encountered: 

b,  cX|+b] 

I(2)  s  j dX[  Jdx2w2(xj,x2) ,  (A-l) 

*2 


where,  specifically  (as  in  section  7.2  of  reference  19), 


(A-2a) 


in  which  x,,x2  are  the  normalized  (i.e.,  "standardized")  random  variables 

X,-X,  X2-Xa  _  -  ' 

Xi  -  — ■ ! - ;  x2  =  —* - ;  p,  2  =  x,x2  =  covan ance  of  x  ; 

°1  CT2 
a2  =  var  Xj  =  (X,  -  Xi)2 ,  etc. 


(A*2b) 


An  analytic  reduction  of  equation  (A-l)  is  sought  from  which  numerical  results  can 
be  obtained,  using  tabulated  functions,  as  well  as  various  analytic  limits.  It  appears  that  in 
the  present  general  form  of  equation  (A-l),  the  best,  i.e.,  formally  simplest,  results  that  can 
be  obtained  involve  at  most  a  single  series  of  tabulated  functions.  Begin  with 

cx,  +bj 

Jw2(xpx2)dx2 

»2 


A-l 


*i  (VAc,  -B/2  VAb  j  )WAb2 

e  ‘  1  r  . 

- /■ . : . e  dz; 

2^1 -P?  2  -Bija^WA., 


p-*f/2  1 

fw  (x  x  )  dx  —  • _ ■ 

0 

xi(ci  ~Pi  2)  +  b2 

+  0 

XlPl  2  a2 

J  W2\X1’X2/UX2  J2n  2 

*2 

V20-p?=) 

.V2(‘-pm1 

(A- 3a) 


where 


(A-3b) 


is  the  familiar  error  integral  in  section  7.1  of  reference  19;  for  associated  tables,  see 
references  20  and  21. 

The  next  step  is  to  consider  the  evaluation  of  integrals  of  the  form 

b,  b,  a, 

I'  =  J dx,  e"x'^0(x,A,  +  Bj  =  JdXj  e"*'^©(x,A1  +  B,)  -  JdXj  e“x'y’20{xIAj  +  B,) ,  (A-4) 

a,  0  0 

which  reduces  to  the  consideration  of  the  generic  relation 


J  s  J  e  ^©(Az  +  §)  dz  , 
0 _ 


(A-5a) 


or  its  equivalent 

eB?/2AJ 

J  —  — — * - 

A 


aA+B  r 

|  e,p|- 


Jt+?y 

2A  T 


0(y) dy , 


(A-5b) 


A  A 

on  change  of  variable  Az  +  B  =  y .  Thus,  using  J  (from  equation  (A-5a))  in  equation  (A-4) 
(from  (A-3))  for  I(2)  in  equation  (A-l): 


I<2)  = 


737 


+I 

2 


b,  a, 

J  ( A,x  +  B, )  e'x,^dx  -  J  ( AjX  +  B, )  ex  J^dx 

0 

j  0(a;x  +  B|)  ex,/2dx  -J  0(  A(x  +  B[)  ex,/2d: 


(A-6a) 


A-2 


with 


a.  =  "fi-ft  j  .  b.  =  -i — bj . ;  a:  =  |  Pi  j  ,  b:  =  ,  •  ,  (A-») 

^O-pH)  V2(‘-pM  V2(>- p-)  '  V2(*-pM 

or,  using  equation  (A-5a),  the  equivalent  is  obtained: 

•••I'2’  =  jjL.{j(b1|A„B,)-J(a,|A„B,)+j(b,|A;,B;)- j(a,|A;,B;)}  .  (A-6c) 

Accordingly,  the  task  here  is  to  evaluate  the  integral  J  (from  equation  A-5a).  Start 
with  the  expansion  of  ©|Az  +  6)  about  B  ,  viz.: 


e(  Az + b) = s(b) + £  e‘(B)i— ^ 

Here,  with  dB  =  dx/V2  and  B  =  x/V2  , 


(A-7a) 


G"(S)  =  )  =  2‘^  =  BV2)  .  n  !.  (A-7b) 

in  which  4>m(x)  =  (l/V2jc)(d,n/dxm)e",,/*  is  tabulated  (see  reference  22),  as  a  form  of  error 
function  and  its  m to-order  derivatives,  here  expressed  as  a  normalized  Gauss  PDF  w,  (x)  = 

,  m  =  0  .  Applying  equation  (A-7)  to  equation  (A-5a)  results  in 


(a|A,B)  =  le-pje(B)  +  |M(-^Jv-'>(V2B) 


(A-8a) 


which  becomes  finally9 


;  .  (Mb) 


9  It  is  shown  in  reference  13  that  the  series,  based  on  equation  (A-7a),  etc.,  converges  only  for  |A|  <  1/V2 
A  complete  treatment  of  (|p)  S  l/V2)  is  provided  in  reference  13,  along  with  numerical  results. 


J 

Ic(fJ;v)  s  — —  f  yv_,e‘y  dy ,  v  =  0  ;  incomplete  J"  -  function  , 
o 


(A-9a) 


This  latter  is  tabulated  also  (see  reference  20).  Here  use 


|e-"*dz=±Jfe(-^); 


Hr  *  ¥ 


Some  limiting  values  of  J^ajA,B  j  are  readily  found: 

j(±Ha.  b) = ±1|s(b) + V2 r(i±>  ><->.(V2b)  .  |a|  £  ^  , 


(A-9b) 


(A-lOa) 


j(±a|oo,B)  =  jV^dz  =  =  j(±a|A,oo) , 

from  equation  (A-5a)  directly,  since  ©(«»)  =  1 : 
j(±°o|oo,B  ,  or  A,«» ,  or  oo,oo  Hf. 

As  a  check,  let  b,,bj  =  °°  and  aj,a2  =  -«> .  1=1,  as  required,  since  w2  is  indeed 

mm 

normalized,  so  that  J  J  w2dx,dy2  =  1 .  Also, 
j(o|A,B)  =  0; 


(A-  10b) 


(A- 10c) 


(A-lla) 


(A-llb) 


J(a|0,0)  =  0; 


(A- 11c) 


and 


where  V2a>(0)  =  ((— 1)“ (2n) !/2B n  !^-n/2tc  (see  equation  (A.1.23)  in  reference  19), 

/.  j(~|a,o)  =  ^^tan'^V2Aj  =  J e'^O^Azj  dz . 

In  addition,  a  further  check  on  the  accuracy  of  the  general  result  (equation  (A-6c)  is 
obtained  by  considering  the  following  special  case: 


Ij2)  *  JjwJ(x1,x2)dx,dx2 


m. 

J 

-  -p  - 

_  j 

0- 

P 

wry  J 

3 

W7j’  J 

+J 

r  - 

l2(u  p!)’ . 

-j 

0 

m 

_£ - r,0 

W?)  J| 

1  1  -i 

~— + — tan 
4  2jc 


Vo-p3) 


=  j_ 
4< 


,  2 
1  +— tan 


n  U(i-p2)\ 


where  J(0| ,  )  =  0  ,  and  where  equations  (A-lOc)  and  (A-l  le)  have  been  used.  Limiting 
cases  of  equation  (A- 12)  here  are 

p  =  l:  w2(x„x2)  =  w1(x1)5(x2-x1)  I,  =  ^  , 
and 

p  =  0:  w2(x„x2)  =  w1(x1)w,(x2)  =  =  7  ’ 

as  expected.  When  0  <  jpj  <  1 ,  then  Ij  obeys  1/4  <  I,  <  1/2  . 

There  remains  the  case  where  A  is  large  (>  3)  but  not  such  that  zA  »  3  ,  z  ~  0 


(A-l  Id) 


(A- lie) 


(A- 12a) 


(A- 12b) 


(A- 12c) 
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Then  the  integeral  for  J  is  split  into  two  parts,  viz.: 


J  = 


’*•  •(**«) 


J+  J  p_I,/20^Az  +  Bjdz  ,  Az0  >  3  /.  9{Az0  j  - 1  -  C^IO-6) ; 

o  x»  J 

=  J  e_zJ^©(  Az  +  B  jdz  + 1  e_I’^dz 


a-  b) + ~  K  vfA )]  • 

a  2:2 

-o  (see  equation  (A  -  9b) ; 

a,(x 

A.ij  ,  since  0^-yLjc  HT6  ,  etc. 

(A- 13) 


(A- 14) 


Accordingly,  with  equation  (A-8b)  applied  to  equation  (A-6b),  (comparatively) 
simple  analytic  expressions  have  been  obtained  for  the  general  probability  integral  I  (see 
equation  (A-l)),  in  terms  of  the  extensively  tabulated  functions  0(x)  in  equation  (A-3a), 
<|>(m)(x)  in  equation  (A-7b),  et  seq.,  and  Ie  in  equation  (A-9a),  namely,  the  error  functions, 
related  derivatives,  and  the  incomplete  F-function.  Computer  evaluation  of  J  and  I  is  then 
straightforward. 
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APPENDIX  B 

CORRECTIONS  FOR  SOME  RESULTS  IN  REFERENCE  3 


In  a  recent  work  (reference  3),  which  among  other  things  considered  the  special 
ternary  case  (M  =  2)  with  purely  (correlated)  Gaussian  noise  and  completely  known  signals 
Si,  S2  (*Si).  i.e.,  coherent  detection,  at  all  signal  levels,  the  limits  on  the  integral 
expressions  for  the  various  decision  probabilities  were  inadvertently  incompletely 
specified.  The  qualitative  nature  of  the  results  were  not  affected,  but  the  quantitative 
expressions  were.  The  correct  analytic  forms  are  given  here. 


Begin  with  equations  (A.l-14a)  and  (A.l-14c),  respectively,  of  reference  3: 

Add  in  the  limits  -  A2  to  y,  +  A,  ;  -A,  to  y2  +  A2 ;} 

Add  in  the  limits  -  A;  to  y2  +  A2 ;  -  A2  to  y,  +  A,  ’  1) 

Also,  for  equation  (A.l-16b), 

-  Yi 

P(D*  *  Pi  I  dyt  |  dy2P?}{y) .  (B-2a) 

-A|  -* 

and,  for  equation  (A.l-17b), 

”  Yj+Aj-A/ 

P'd^Pj  /dy,  Jdy,Pf(y).  (B-2b) 

-a2 


The  results  of  these  limit  modifications  are  that  equation  (A.l-33b)  becomes  (with  d>  a  sk'ls) : 


pg'-=a-[i+e(B®)], 

' 

B(i>  a  max 

mu 

aGV<>  .  logit,  a0Vd> 

logW 

UJ 

Itt +V2a0v*  •  ~ir 

V2a0Vd> 

(B-3) 


while  equation  (A.l-33c)  is  replaced  by  equation  (A.l-14a)  (see  equations  B-l,  etc.). 


B-l 


Similarly,  for  equation  (A.l-34b): 


p(Dl'=JHi+e(B£,)],with 


b!2. 


max 


a0V<fr  ,  log(i2  a0Vfc 

2^2  V2acV$  *  V2an^ 


logfe) 


(B-4) 


In  the  same  way,  for  equation  (A.  1-35): 


R*  - 

M“2_  2 

l-q@ 

r  min  > 

r  k 

2V2  V2r 

-p^b'D-p^bSS,) 

,  T  sa0V<b  ; 

(B-5) 

> 

so  that,  for  equation  (A.  1-36), 


r-4,.*3  ■ 

which  gives,  for  equation  (A.  1-27), 

AR“-!=^H^)+20(:72)]' 

Finally,  equation  (A.  1-38)  becomes,  for  q  =  p,  =  p2  =  1/3 


(B-6) 


(B-7) 


(B-8) 
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